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Abstract 

Thin film junctions and bilayers of the doped topological insulator Bi2Se3 and the s-wave super- 
conductor NbN were found to exhibit conductance spectra with prominent zero bias and coherence 
peaks. Various tunneling models with different pair potentials have failed to fit our data, except 
for the triplet px + ipy pair potential, which breaks time reversal symmetry, but yielded reasonably 
good fits. This provides evidence for proximity induced triplet superconductivity in the Bi2Se3 
layer near the interface with the NbN film. 
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Topological superconductors (TSC) are interesting since they are predicted to support 
Majorana fermions (MF) which are protected against disorder and decoherence, and are 
supposed to play a role in future quantum computers [lH3]. This type of superconductors 
can be realized by either doping of a topological insulator (TI) [IJ-E], or by putting a TI 
or a semiconductor with strong spin orbit coupling in good contact with a conventional 
superconductor In the former case, an intrinsic bulk superconductivity occurs, while 

in the second case, superconductivity is induced by the proximity effect. Recently, sup- 
porting experimental evidence for the existence of TSC and MF bound states was found in 
superconducting nano- wires of the proximity induced hybrids of InSb/NbTiN and InAs/Al 
[Tin ITT] - In junctions of other TSC materials such as CuxBi2Se3, Bi2Se3 — Sn, Bi2Te3 — Bi 
and Bi2Se3 — NbN with normal metals, generally robust zero bias conductance peaks 
(ZBCP) were observed, but conclusive evidence that these are due to MF in a TSC rather 
than Andreev bound states (ABS) in an unconventional superconductor, is still lacking 
[SHH]. In the present study we extended our previous investigation of Au — Bi2Se3 — NbN 
junctions of large overlap area (~ 100 x 150 yum^), to smaller ramp-type junctions of a 
much smaller area (5 x 0.5 /xm^), and to even smaller area junctions of a few nm diameter, 
resulting from crashing of an STM tip into Bi2Se3 — NbN bilayers. The observed ZBCP 
disappeared above a proximity induced transition to superconductivity at Tc(junction)~8 
K, while the coherence peaks survived up to Tc of the NbN electrode at ~10 K. The 
conductance spectra at low temperature, which showed a ZBCP and coherence peaks, could 
be fitted very well with a minimal number of fitting parameters only when using the triplet 
and spinfull + ipy pair potential, which breaks time reversal symmetry. This unexpected 
result calls for further investigation into the spin state nature of the pair potential at the 
Bi2Se3 — NbN interface. We note that indications for spin-triplet superconductivity were 
found recently also in single crystals of the doped topological insulator CuxBi2Se3 in high 
pressure and point contact measurements [121 US] • 

Preparation and characterization of the NbN, Au and Bi2Se3 thin films by laser ablation 
deposition was described in our previous study [9]. All films and bilayers were deposited 
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on (100) SrTiOa (STO) wafers of 10 x 10 mm^ area. Here we shall only elaborate on the 
preparation of the fully photographically patterned ramp-type junctions, whose schematic 
cross-section is shown in the inset of Fig. 4, and are also described in |14j. First, a bilayer 
of the base electrode consisting of an insulating and smooth STO layer of 50 nm thickness 
was deposited under vacuum at 150° C on a 70 nm thick NbN layer. This bilayer was 
then patterned in the form of Fig. 2(c) of Ref. [2] by Ar ion milling on half the wafer, 
and left in ambient air for a day to create the 1-2 nm thick tunneling barrier of the 
native oxide layer ("nox" - marked by the inclined red line in the inset of Fig. 4) [T5] . 
Then the cover electrode of 100 nm Au on 70 nm Bi2Se3 was deposited and patterned 
to produce ten ramp junctions on the wafer, each with about 5 x 0.5 /xm^ junction area. 
Note that no current could flow via the thick STO layer in the 5x2 yum^ overlap area of 
the cover and base electrodes in the direction normal to the wafer. Conductance spectra 
were measured by the four-probe technique simultaneously on all the ten junctions, using a 
4x 10 array of spring loaded gold coated tips pressed on the 40 gold contact pads on the wafer. 

First, we demonstrate that ZBCPs in the junctions were observed only below the 
transition temperature of the proximity induced superconducting layer of the Bi2Se3 in 
the interface region with the NbN film. Fig. 1 shows normalized conductance spectra at 
different temperatures of a large junction with overlap area of 100 x 30 fim^, together with 
the resistance versus temperature of this junction in the inset. Below the transition to 
superconductivity of the NbN electrode at Tc{NbN) ~10 K, one can see a clear second 
transition at about 8 K which is due to the proximity effect in the Bi2Se3 layer near the 
interface with the NbN layer. This second transition therefore originates in the junction 
and is denoted by Tc{junction). At low temperature the ZBCP is robust (not shown). 
With increasing temperature it decreases in height but persists up to Tc{junction) ~8 K, 
above which it disappears and only the weak coherence peaks remain up to of the NbN 
electrode. A clear correlation is therefore shown between the presence of ZBCPs in the 
main panel and the proximity induced superconductivity in the junction. 
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FIG. 1: (Color online) Normalized conductance spectra of a large junction with 100 x 30 ;um^ 
overlap area at various temperatures. The overlapping layers of this junction (as in Fig. 2 (a) of 
Ref. [9j) are: 100 nm Au on 20 nm Bi2Se3 on 1-2 nm native oxide barrier on 70 nm NbN. The inset 
shows the resistance versus temperature of this junction, with the transition of the NbN electrode 
at Tc{NbN) and the proximity induced transition of the junction at Tc{junction). 

Next, we use a conductance spectrum with a robust ZBCP and coherence peaks from 
our previous study (Fig. 5(d) of Ref. [^) as a test case for finding which pair potential 
(PP) fits this data best with a minimum number of fitting parameters. Fitting was done 
using the standard Blonder, Tinkham and Klapwijk (BTK) model [16], with modifications 
of it to other PP symmetries as needed [T7HT9] . The standard parameters of the fits were 
the barrier strength Z, the pairing amplitude or energy gap Aq, and the finite lifetime 
broadening F. Various pair potentials were used, and the results of the fits are shown in 
Fig. 2 together with the original experimental data. Singlet s-wave pair potential could fit 
the coherence peaks but not the ZBCP, while a fully gapped TSC, as well as a TSC with 
point nodes, could barely fit the ZBCP and the two dips beside it, but not the coherence 
peaks [19] . We note that since there is no direct contact between the NbN electrode and the 
gold overlayer, the observed coherence peaks in Fig. 2 could not originate in any NbN-Au 
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FIG. 2: (Color online) Normalized conductance spectra of a large junction of a previous study 
of our group (Fig. 5(d) of Ref. [9j), together with fits to different theoretical models and pair 
potentials. 



junction. It is also noted that the conductance spectra predicted by Yamakage et al. [19] 
who used various PP as proposed by Fu and Berg [20], had never shown coherence peaks 
besides the ZBCP. Moreover, if we assume the existence of two bands of the surface and 
bulk states in the superconducting Bi2Se3 near the interface, with two different and 
A;, gaps, we could fit the data of Fig. 2 quite well using the s-wave PP. This however 
necessitates very different barrier strengths (Z values) for the two bands, which is unphysical 
and thus rules out the s-wave fits. We therefor tried the singlet d-wave PPs d^2_y2 and 
dx2-y2 + idxy [17|, ignoring the hexagonal symmetry of Bi2Se3, but compensating for it 
by taking weighted tunneling contributions from the node and antinode directions while 
keeping the same fitting parameters Z, Aq and F for both orientations. This weighted 
sum of the calculated individual conductances along the anti-node G(0°) and node G(45°) 
directions amounts to a[G{0^) — 1] + b[G{A5^) — 1] + 1 where a and b are fit parameters. 
They represent the relative weight of the two orientations b/a, and the normalization factor 
to the experimental data a. Fig. 2 shows that the fits improved, but not sufficiently to 
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reproduce the whole width of the measured ZBCP. In passing we add that ignoring the 
hexagonal symmetry of the Bi2Se3 layer on the cubic NbN film can be further justified 
by noting its mosaic structure as shown in Fig. 12 of Ref. [9]. Finally, we tried the odd 
parity triplet PP + ipy which breaks time reversal symmetry, but shows consistency with 
some experimental results [121 [13]. Following [18], we used Ai-|- = AQsin6{cos(f) + sincp) and 
A-f-f- = AQsin6{cos(f) + isincj)) with 9 and being the polar and azimuthal angles, for the 
two (tetragonal) symmetries Eu(l) and Eu(2), respectively, while Aj^^ = A-f-^^ = A^^-f- = 0. As 
before, we sum over weighted contributions to the conductance from two different interfaces 
(xy and yz, where xy represents the a-b plane of the Bi2Se3 layers), and keep the same 
basic fitting parameters Z, Aq and P for both interfaces. The justification for taking the 
weighted contributions of the two interfaces is that both are present in our junctions, either 
via the roughness of the interfaces or due too the inclined geometry of the barrier in the 
ramp junctions. Fig. 2 shows that while the best fit using the -E'u(l) PP fails to fit the 
ZBCP, the Eu(2) PP fits the data quite well. All the features of the spectrum are captured 
satisfactorily and only small deviations are found which can be attributed to the tetragonal 
rather than hexagonal symmetry used in the model calculations. It is thus concluded that 
if this good fit using the triplet px + ipy wave PP has any physical significance rather than 
just being coincidental, it should also fit conductance spectra of other Bi2Se3-NbN junctions 
of various types. In the following we show that this is actually the case in both ramp type 
junctions and point contact junctions of this system. 

Figs. 3 and 4 show normalized conductance spectra obtained on ramp junctions whose 
cross section is shown schematically in the inset of Fig. 4. At low temperature and zero 
magnetic field. Fig. 3 shows tunneling-like spectra with robust ZBCP and coherence peaks 
similar to those of Fig. 2, while Fig. 4 shows Andreev-like spectra, but still with clear 
ZBCP and coherence peaks. Both spectra could be fitted reasonably well using the triplet 
EuC^) Px + ipy PP as before, but with different barrier strengths values Z (1.5 and 0.55) in 
accordance with the tunneling-like and Andreev-like behaviors. Fig. 3 also shows spectra 
at different temperatures which show interesting correlation with the resistance versus 
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FIG. 3: (Color online) Normalized conductance spectra of a small ramp type junction of 5 x 0.5 iJ,m? 
area at various temperatures, together with a p-wave fit to the data at 1.86 K. The inset shows 
the resistance versus temperature of this junctions, with the kinks at Td and Tc2 indicating the 
temperatures at which the narrow and broad ZBCPs disappear, respectively. 
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FIG. 4: (Color online) Normalized conductance spectra at 1.82 K and different magnetic fields 
of another ramp junction on the same wafer as that of Fig. 3. This junction is more transparent 
than that of Fig. 3, as the p-wave fit to the data at 1.82 K shows. The inset depicts a schematic 
cross section of a typical ramp junction. 
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FIG. 5: (Color online) Two normalized conductance spectra at 4.2 K of a point contact junction 
of a Pt/10%Ir STM tip which was crashed into a bilayer of 10 nm Bi2Se3 on top of 70 nm NbN. 
The curve is a fit to the Eu{l) p-wave model. 

temperature results of the inset to this figure where two kinks are found at Td ~8 K and 
Tc2 ~10 K. The robust ZBCP of ~1 mV width persists up to about T^i, while above it 
and up to Tc2 only a broader ZBCP of ~2 mV width remains. This is apparently due 
to the merging of the decaying ZBCP and coherence peaks with increasing temperature. 
Fig. 4 shows that with increasing magnetic field at low temperature the spectra become 
narrower, decay in height, and smeared as the ZBCP and coherence peaks merge with one 
another. In all our junctions we had never seen a ZBCP increasing with field as in Refs. 
pm [TT] . We note that while the Eu{l) PP failed to fit the tunneling spectra of Fig. 3 
(similarly to the case of Fig. 2), it did fit the data of Fig. 4 equally well as the Eu(2) PP 
with the same Aq and F values, but with a lower Z value (0.2 versus 0.55, not shown). 
Another observation is that the gap energy of the Andreev-like spectra of Fig. 4 (2.1 mV) 
is quite larger than that of the tunneling-like spectra of Fig. 3 (1.3 mV). This result seems 
to be due to the larger drop in the gap energy at the interface in the less transparent 
junction. From the good fits of the data in Figs. 3 and 4 we conclude that also in small 
ramp junctions with different transparencies, the triplet Px+ipy PP describes our results well. 
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Fig. 5 depicts two normalized conductance spectra at 4.2 K of a point contact junction 
on an in-situ prepared bilayer of a 10 nm thick Bi2Se3 layer on top of a 70 nm NbN film 
(no native oxide at the Bi2Se3-NbN interface). The wafer consisted of two halves, one 
half with the bilayer and the other half with a reference bare film of 70 nm NbN. In the 
STM mode (vacuum gap between tip and sample), the tunneling spectra measured on both 
halves did not reveal clear signatures of superconductivity, apparently due to a high Z oxide 
layer on the NbN, and the presence of a highly deteriorated top layer of the Bi2Se3 which 
developed when the sample was exposed to ambient air. Only when the Pt/10%Ir STM tip 
was intentionally crashed into the top of the Bi2Se3 layer of the bilayer, forming a point 
contact of a few nm in diameter, the spectra of Fig. 5 could be measured. Here again use 
of the triplet px + iPy PP fitted the data very well, but now with the -£'«(!) PP and a Z=5 
value, the highest used in the present study. Attempt to fit the data of Fig. 5 also with the 
Eu(2) PP using the same Z, Aq and F parameters failed, as the three main features of the 
spectrum (the ZBCP maximum, dips minima, and coherence peaks maxima) could not be 
fitted all at once (not shown). Here again, we come to the same conclusion that the triplet 
Px + "iPy PP describes our data well, also in this very small point contact junction of a few 
nm diameter. 

Fu and Kane had discussed the issue of a proximity effect between an s-wave supercon- 
ductor and a topological insulator pT|. They predicted that the superconductor induces in 
the surface states of the TI a superconducting PP which resembles that of a spinless px + ipy 
state. By a proper transformation of the base functions they showed that this PP could be 
made precisely spinless, and with it the Hamiltonian preserves time reversal symmetry. The 
presently used triplet Px + ipy PP which fits our data nicely is clearly a spinfull PP which 
breaks time reversal symmetry. To reconcile these two seemingly contradictory results, we 
propose a scenario in which the odd-parity, interorbital spin-triplet A4 PP of [20j, which 
breaks time reversal symmetry, is induced in the Bi2Se3 near the interface, rather than 
the A2 PP which preserves it and is supposed to exist in CuxBi2Se3. A4 is a spinfull PP 
which is consistent with our present results. It is possible that the proximity induced A4 
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PP originates in the bulk states at the chemical potential of our highly doped Bi2Se3 films, 
and if its contribution to the tunneling conductance at the interface is dominant, it might 
mask the surface contribution that could come from A2. We note that in the proximity 
region of the interface both A4 and A2 can coexist, while in a TSC such as possibly 
CuxBi2Se3, they could not since they belong to different irreducible representations. The 
observed ZBCPs are thus due to a spinfuU triplet PP and not to Majorana bound states as 
previously envisioned. Wc believe that if the bulk conductivity of our Bi2Se3 films could 
be significantly reduced, observation of MF in our junctions could be realized. We also 
note that in most worldwide films and single crystals of Bi2Se3 the bulk doping with Se 
vacancies is a major problem, which yields dominant metallic bulk conductance that masks 
the topological surface conductance. The later of course, could yield MF bound states 
when in the superconducting phase. In the near future, we plan to use low energy muon 
scattering from bilayers of Bi2Se3 on NbN in order to see if there is a net magnetization 
at zero field in the interface layer, originating in the same-spin triplet PP observed in the 
present study. 

To summarize, we found a signature of proximity induced triplet superconductivity by 
NbN in the doped topological insulator Bi2Sc3. Conductance spectra of various junctions 
with different transparencies that have robust ZBCPs and coherence peaks, could be fitted 
quite well using the time reversal symmetry breaking pair potential + ipy Thus, the 
observed ZBCPs of this odd parity, same-spin triplet pairs represent zero energy surface 
bound states which do not originate in Majorana fermions. 
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